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Abstract 

The first-order irreversible phase transitions (IPT) of the Yaldran-Khan model 
(Yaldran-Khan, J. Catal. 131, 369, 1991) for the CO + NO reaction is studied 
using the constant coverage (CC) ensemble and performing epidemic simulations. 
The CC method allows the study of hysteretic effects close to coexistence as well as 
the location of both the upper spinodal point and the coexistence point. Epidemic 
studies show that at coexistence the number of active sites decreases according to a 
(short-time) power law followed by a (long-time) exponential decay. It is concluded 
that first-order IPT's share many characteristic of their reversible counterparts, such 
as the development of short ranged correlations, hysteretic effects, metastabilities, 
etc. 

1 Introduction 

The study of critical phenomena occurring in adsorbed overlayers is a topic that has 
attracted increasing attention [1]. In fact, the understanding of the behavior of atoms 
and molecules absorbed on different surfaces is essential for many branches of science 
(chemistry, physics, biology, etc). Within this wide context, surface phenomena occurring 
under non-equilibrium conditions are far from being understood. Therefore, they have 
become the center of great attention [2]. In particular, irreversible catalytic reactions 
exhibit a very rich and complex behavior that includes oscillations, bifurcations, chaos 
metastability, irreversible phase transitions (IPT's), etc. [3, 4, 5, 6, 7, 8, 9]. 

The aim of this work is to perform an extensive numerical study of the first-order IPT 
characteristic of a model for the NO + CO reaction early proposed by Yaldran and Khan 
(YK) [4]. IPT's in reactive systems take place between an active regime with sustained 
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outcome of the reaction product from the catalytic surface and an absorbing (or poisoned) 
state where the catalyst becomes fully covered by one or more types of reactants. Since 
the systems cannot escape from the absorbing state the transitions are irreversible [2, 5, 6]. 
The study of IPX's has gained increasing attention in the field of nonequilibrium statistical 
physics since the pioneering work of Ziff, Gulari and Barshad (ZGB) that introduced a 
simple lattice model for the catalytic oxidation of CO [7]. Particularly interesting arc first- 
order IPX's that can be characterized by an abrupt change in the density of reactants and 
in the rate of production due to a tiny change in the (external) control parameter, which 
is usually the pressure. In fact, evidence of such kind of transition has been reported 
in catalyzed reactions [8, 9]. Xhe occurrence of hysteretic effects around a first-order 
transition point has also been observed [10]. Very recently, we have performed a detailed 
study of the first-order IPX of the ZGB model presenting conclusive evidence of hysteretic 
effects [11]. As anticipated above, the aim of this work is to characterize the first order 
IPX of the YK model [4] for the CO + NO reactions by means of Monte Carlo simulations 
with particular emphasis on the occurrence of hysteretic effects and by performing useful 
comparisons with our previous study of the ZGB model [11]. 

Xhe manuscript is organized as follows: in Section II we describe the YK model and 
three different simulation methods used in the work, namely the standard ensemble, the 
constant coverage ensemble and the epidemic approach. In Section III we present and 
discuss the obtained results, while our conclusions are stated in Section IV. 



2 The model and Monte Carlo simulation methods. 

2.1 The YK model for the CO + NO reaction. 



Yaldran and Khan [4] have proposed a lattice gas model for the catalytic reaction of 
CO + NO based on the Langmuir-Hinshelwood mechanism. Xhe reaction steps are as 
follows: 



where S represents an unoccupied site on the catalyst surface, 2S represents a nearest 
neighbor (NN) pair of such sites, (g) indicates a molecule in the gas phase and (a) indicates 
an species adsorbed on the catalyst. Xhe reactions given by equations (3) and (4) are 
assumed to be instantaneous (infinity reaction rate limit) while the limiting steps are 
the adsorption events given by equations (1) and (2). Xhe YK model is similar to the 
ZGB model for the CO + O2 reaction, except that the O2 is replaced by NO, and NN 



NO{g) + 2S^ N{a) + 0{a) 

CO{g) + S ^ CO(a) 
CO{a) + 0(a) ^ C02ig) + 2S 
N{a) + N{a) ^ N2{g) + 2S 



(1) 
(2) 
(3) 
(4) 
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atoms, as well as NN CO — O pairs, react. For further details on the YK model see 
[4, 12, 13, 14, 15]. 

2.2 Monte Carlo simulation method using the standard ensem- 
ble 

Monte Carlo simulations are performed on the hexagonal (triangular) lattice of side L, 
assuming periodic boundary conditions. In the standard ensemble the procedure is as 
follow: let Pno and Pco be the relative impingement rates for NO and CO, respectively, 
which are taken to be proportional to their partial pressures in the gas phase. Taking 
Pco + Pno = li such normalization implies that the YK model has a single parameter 
that is usually taken to be Pco- CO and NO adsorption events are selected at random 
with probabilities Pco and 1 — Pco, respectively. Subsequently, an empty site of the 
lattice is also selected at random. If the selected species is CO, the adsorption on the 
empty site occurs according to equation (2). If the selected molecule is NO, a NN site 
of the previously selected one is also chosen at random, and if such site is empty the 
adsorption event takes place according to equation (1). Of course, if the NN chosen site 
is occupied the adsorption trial is rejected. After each successful adsorption event all NN 
sites of the adsorbed species are checked at random for the occurrence of the reaction 
events described by equations (3) and (4). 

The Monte Carlo time step (MCS) involves adsorption attempts, so that every site 
of the lattice is selected once, on average. Simulations are started with empty lattices. 
The first 10^ MCS are disregarded to allow the system to reach the stationary regime and 
subsequently, averages are taken over 4 x 10^ MCS. During the simulations, the coverages 
with CO, O and {Oco, do and On, respectively) as well as the rate of production of CO2 
and N2 {RcO'zi Pn2, respectively) are measured. The phase diagram of the YK model is 
similar to that of the ZGB model [7], in the sense that both second- and first- order IPT's 
are observed. However, in contrast to the ZGB model where the absorbing (poisoned) 
states are unique, in the case of the YK such states are mixtures of 0(a) -|- N{a) and 
CO{a) -\- N{a) as follows from the observation of the left and right sides of the phase 
diagram, respectively (figure 1(a)). 

The IPT observed close to P^,q = 0.184(1) is continuous and therefore of second-order 
(see figure 1). Our estimation of is in agreement with previous calculations, namely 
P^o = 0.185(5) (reference [12]) and P^q = 0.185(2) (reference[4]). 

More interesting, an abrupt first-order IPT is also observed close to P^q = 0.3545(5) 
(figure 1(a) and (b)), in agreement with previous simulations (references [4, 12]). It should 
be noticed that a reactive window for the YK model is observed on the hexagonal lattice 
(see figure 1) while such window is absent on the square lattice [4, 12], pointing out the 
relevance of the coordination number on the reactivity. 
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Figure 1: Phase diagram of the YK model on the hexagonal lattice of size L=128 LU. (a) 
Plots of ^c7o(A), ^o(n) and ^jv(O) versus < Pco >■ (b) Plots of RNii^) and i^coalO); 
measured in units of number of and CO2 molecules removed from the lattice per unit 
of area and time, respectively; versus < Pco >■ 

2.3 The constant coverage ensemble 

The constant coverage ensemble (CC) was early introduced by Ziff and Brosilow [16] 
to study the first-order IPT of the ZGB model. The first step of the CC ensemble is 
to obtain a stationary configuration using the standard ensemble. Then the system is 
actually switched to the CC ensemble, where the density 6co is now kept as constant 
as possible. For this purpose, if Oqo is below the established value, CO is adsorbed 
on a randomly selected site. Otherwise, if 9co is greater than the desired value, NO 
adsorption attempts on randomly selected sites are performed. Now the effective CO- 
pressure (< Pco >) is given by the ratio of CO-adsorption attempts to the total number 
of adsorption attempts. 

It is worth mentioning that using the CC method as proposed by Ziff and Brosilow 
[16], the actual coverage with CO is not strictly constant but it is affected by fluctuations 
of the order of fa Very recently [17], a true constant coverage ensemble, where the 

number of particles remains strictly constant, has been proposed and used to study the 
conserved contact process. This new method seems to be very useful for the study of 
single particle systems, such as the contact process, branching annihilating walkers, etc., 
although their implementation for complex multiparticle systems, such as the YK model, 
does not appear to be straightforward. 
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In simulations, < Pco > is averaged over tm time steps. Subsequently, 6co is in- 
creased stepwise and measurements of < Pco > are performed after some waiting time 
Tw to allow the relaxation of the system. In this way the growing branch (GB) of the CC 
loop can be obtained. Then, after reaching a large value of 9co {^co ~ 0.88 in this work 
in order to prevent the system from reaching an absorbing state as shown in figure 1), 
dco is decreased stepwise. Using this procedure the decreasing branch (DB) of the CC 
loop can be recorded. It should be noticed that in the CC ensemble 6co assumes the role 
of the control parameter. Further details on the CC method can be found in references 
[6, 11]. 

2.4 Epidemic simulations close to coexistence. 

Another powerful approach to the study of IPX's is to perform the so-called epidemic 

analysis (EA) [18, 19]. In EA the simulations start from a configuration very close to the 
absorbing state and subsequently the time evolution of the system is followed using the 
standard ensemble (SE). In order to generate the initial configuration, a natural absorbing 
state has to be achieved first using the SE and following the dynamics of the system. This 
procedure assures the development of the characteristic correlations among the reactants. 
Taking the obtained absorbing state, a small patch of empty sites is created close to the 
center of the sample. Subsequently, during the time evolution of the system the average 
number of empty sites {N{t)) and the survival probability P{t) of the active state are 
recorded. Each single EA stops when the system is trapped in the absorbing state so that 
N{t) — 0. In order to obtain reliable results, the quantities of interest have to be averaged 
over a large number of independent EA (~3xl0^ runs in the present work). 

Performing EA close the second-order IPT the scaling ansatz N{t) oc t'', where r] 
is an exponent, has been proposed to hold at criticality [18, 19]. This observation is 
in agreement with well established ideas developed by studying equilibrium (reversible) 
phase transition: scale invariance reflects the existence of a divergency in the correlation 
length. However, close to first-order transition it is also well known that correlations are 
short ranged, preventing the observation of scale invariance. Recently Monetti and Albano 
[11] have proposed that, at coexistence, N{t) should decrease according to a short-time 
power law followed by a long-time exponential decay, so that: 

A^Woc(^)-''^^^exp[-(l)] (5) 

where T sets a characteristic crossover time scale and r^^ff is an effective exponent. It 
should be noticed that equation (5) holds for the ZGB model at coexistence [11]. 
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Figure 2: Plots of 9co versus < Pco > obtained using the CC ensemble with L — 32 LU, 
tr — 50 MCS and tm — 50 MCS. The arrows show the narrow vertical region of the loops 
(VR) and the upper spinodal point {Pco)^ respectively. More details in the text. 

3 Results and discussion. 

3.1 Study of the hysteretic effects using the CC ensemble. 

In order to study hysteretic effects CC simulations using lattices of different sizes and 
various values of tw have been performed. Since the measurement in time (tm) has to be 
reduced in order to observe hysteretic effects, results are averaged over several (10 — 100) 
different loops, depending on the lattice size, in order to obtain acceptable statistics. 

Figure 2 shows a plot of 9co vs. < Pco > as obtained using L — 32LU, tw = 50 
MCS, and tm = 50 MCS. For this small lattice size the relaxation time is quite short, so 
that hysteretic effects are absent. This result is in agreement with similar measurements 
performed applying the CC ensemble to the ZGB model [11]. Note that a narrow vertical 
region close to the center of the loop for (< Pqo > ~ 0.35) can also be observed in figure 
2. On the other hand, the L-dependent upper spinodal point {Pco) ^^"^ clearly be 
observed as shown in figure 2. On increasing the lattice size (L = 6ALU, in figure 3), 
the onset of hysteretic effects can be observed for tw = 75 MCS (figure 3(a)), while such 
effects become almost neghgible if the waiting time is increased (figure 3(b) for tw = 400 
MCS). Furthermore, as in the case of figure 2, a vertical region located at the center of 
the loop and slightly above < Pco > ~ 0.35 can be observed in figure 3(a). Also notice 
that such region becomes well defined when tw is increased (figure 3(b)). 

On increasing the lattice the hysteretic effects can be observed even taking larger 
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Figure 3: Plots of 9co versus < Pco > obtained using the CC ensemble with L = 64 
LU, and (a) tw = 75 MCS, tm = 50 MCS and (b) tw = 400 MCS, tm = 100 MCS. The 
arrows pointing up and down, in figure (a), show the growing and decreasing branches 
(GB and DB) of the loop, respectively. More details in the text. 



values of T]^. In fact, for L = 128LU (figure 4(a)) and L = 256 (figure 4(b)), hysteresis 
is still observed for tw = 1700 MCS and tw = 2700 MCS respectively. 

It should be noticed that the loop becomes narrow when the lattice size is increased, 
as follows from the comparison of figures (3) and (4). Comparing these figures it also 
becomes evident that while the location of Pqq is shifted systematically toward lower 
values when L is increased, the location of the vertical region (close to the center of the 
loops) remains almost fixed very close to Pco = 0.3515 (see figure 4(b)). 

The evaluation of a CC loop for L — 1024LC/ requires huge CPU resources, so we have 
restricted ourselves to the case tw — 600 MCS and tm — 700 MCS, since considerably 
larger values of Tm became prohibitive. In this case (figure 5) hysteretic effects are quite 
evident and also, the growing and decreasing branches of the loops arc almost vertical. 
However, the vertical region at the center of the loop, previously observed using smaller 
lattices, is no longer found. 

In order to gain insight into the behavior of the system close to coexistence it is useful 
to analyze snapshot configurations, as shown in figure 6. Snapshots are obtained for some 
relevant points as indicated in figure 4(a). Just at the upper spinodal point (figure 6(a)) 
one observes that the active phase remains homogeneous but the nucleation of few CO- 
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Figure 4: Plots of 9co versus < Pco > obtained using the CC ensemble and taking: (a) 
L = 128 LU, Tw = 1700 MCS, and tm = 300 MCS. The points labeled A, B, C, D and E 
correspond to the coverages used to obtain the snapshot configurations shown in figure 6 
(a)-(e), respectively, (b) L = 256 LU and tw = 2700 MCS and tm = 300 MCS. 
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Figure 5: Plots of 6co versus < Pco > obtained using the CC ensemble for (a) L — 1024 
LU, and tw = 600 MCS, and tm = 700 MCS. 
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clusters has already started. The biggest one (close to the lower right hand side in figure 
6(a)) is compact (except for the presence of few A^- atoms embedded into the bulk of the 
CO-cluster) and has reached the critical nucleation size. At the upper spinodal point one 
has 6co ~ 0.045. On increasing the CO coverage this critical nucleus grows up into a 
solid and compact CO cluster surrounded by the (homogeneous) active phase, as shown 
in figure 6(b). The coverage with CO has increased up to Oqo ~ 0.28, the compact CO 
cluster does not percolate and its interface is essentially convex. Further increasing 9co 
causes the CO cluster to percolate along one direction of the sample, figure 6(c). The 
hysteresis loop also changes for the growing branch (as in figure 6(b)) to the central region 
mih. 9co ~ 0.63 in figure 6(c). The percolating cluster is quite stable and has an interface 
essentially fiat with a length of the order of 2L. An additional growth of CO coverage 
causes the CO cluster to percolate along both directions of the lattice (figure 6(d)). Here 
the interface of the cluster is essentially concave and the active phase is surrounded by 
the massive cluster. Since the coverage with CO is the same {9co ~ 0.63) in both figures 
6(c) and 6(d), the presence of hysteretic effects has to be related to the curvature of the 
interface of the CO cluster. Subsequently, on decreasing 9co one observes the formation 
of a CO cluster that percolates along only one direction of the sample (figure 6(e)). Such 
kind of clusters are observed along the vertical region and are characterized by an almost 
flat interface with an infinite effective curvature radius. Notice that despite the fact that 
the CO coverage is the same in both figures 6(b) and 6(e) {9co ~ 0.28), the curvature of 
the interface of the CO cluster is quite different and hysteretic effects are observed (see 
figure 4(a)). 

Summing up, the sequence of snapshots shown in figure 6, as well as additional figures 
not shown here for the sake of space, allows us to infer the following behavior of the 
interface of the CO cluster along the hysteresis loop: between the spinodal point and 
the growing branch the curvature radius of the interface is positive and finite. Along 
the vertical region the radius of curvature is infinite, while within the decreasing branch 
the radius of curvature is negative and finite. Almost horizontal regions observed in the 
loops correspond to regions where the curvature radius crosses over between two different 
behaviors. 

In order to quantitatively analyze the behavior of Pco corresponding to the different 
branches and the vertical region, it is assumed that the growing (decreasing) branch starts 
at the point where both curves merge (split out). Figure 7 shows the dependence of the 
location of the growing branch and the decreasing branch {Pco Pco ^ respectively) 
on the inverse of the lattice size. The L— dependence of Pco the vertical region {Pco) 
has also been included for the sake of comparison. 

As already discussed the location of all relevant points, namely Pco with X = GB, DB 
and VR depend on the curvature radius (s) of interface of the massive CO cluster. Such 
dependence can be written as follows: 

P^o-Pco{L^^) + F''{s), (6) 
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Figure 6: Typical snapshot configurations obtained using the CC ensemble close to coexis- 
tence for L = 128 LU and the points A, B, C, D and E shown in figure 4(a). (a) Snapshot 
corresponding to the upper spinodal point S=A, (b) and (d) snapshots corresponding to 
the points B and D respectively, (c) and (e) CO-percolating clusters obtained in the 
central region (points C and E respectively). Here o denotes 0(ad), +denotes A^(ad) and 
• denotes CO (ad). 
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Figure 7: Plots of < Pco > versus measured in the growing branch (O); decreasing 
branch (A), and the vertical region (□). The straight lines correspond to the best fits of 
the data that extrapolate to L — > oo. L is measured in LU. 

where Pqq{L oo) is the location of the point under consideration after proper extrap- 
olation to the thermodynamic limit and F{s) is an s-dependent function. For the vertical 
region one has s ^ oo and P^^q is almost independent of L, so F^-^(oo) 0, as shown 
in figure 7. In contrast, for the DB and the GB , s is finite and of the order of — 1/L 
and 1/L, respectively. So, one has F^^{s) ^ —A/L while F'^^{s) ^ B/L. AH these 
arguments can be confirmed by the results shown in figure 7 and the extrapolated points 
are: 



respectively. Also, A ^ 0.215(5) and B,Ki 0.12(2) are obtained. 

The observed behavior allows us to identify Pcoi^ ~^ °°) ^ coexistence point 
pCoe^' ^ 0.35145(5) in excellent agreement with the value Pco = 0.35140(1) reported by 
Brosilow and Ziff ([12]). 

This result is in contrast with measurements performed with the ZGB model. In fact, 
for the ZGB systems the vertical region is not observed while the locations of the growing 
and decreasing branches are almost independent of the lattice size [11]. Consequently, the 
CC ensemble does not provide a method suitable for the location of the coexistence point 
that has to be estimated using the spontaneous creation method [11]. It is expected that 



P^^(L^oo) = 0.3514(3), 
P^f{L oo) = 0.3517(3), 
P^^{L ^ oo) = 0.35145(5), 
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Figure 8: Plots of the local CO— coverage fluctuations {P{9co)) versus 9co measured 
using the constant coverage ensemble keeping the global CO coverage at 9co — 0.040 and 
using lattices of different sizes as indicated in the figure, where L is measured in LU. The 
inset shows a detailed view of the right-hand side of the distributions. Results averaged 
over 3 x 10^ different configurations. 



the difference observed between the models may be due to the different behavior of the 
interface of the massive CO cluster. So, we are planning to perform extensive simulations 
on this subject to clarify this open question. 

Pointing our attention to the upper spinodal point, it is found that its location depends 
on the lattice size, so P^q{L) can be determined from the loop, as shown in figures (2- 
4). It is expected that this dependence of P(jq{L) is due to local fiuctuations of the CO 
coverage that take place during the nucleation of the critical cluster. In order to check this 
conjecture, the probability distribution of the local CO coverage {PiOco)) m small patches 
of side Lo = 30LU was measured for lattices of different sizes. It should be noticed that, 
in order to perform these measurements, the CC algorithm that keeps the global coverage 
of CO almost constant has been used. Setting 9co = 0.040, i.e. a value slightly smaller 
than the coverage at the upper spinodal point, the probability distributions shown in 
figure 8 have been obtained. For a rather small lattice (L = 48 in figure 8) P{Oco) is 
almost symmetric around the maximum 9"^^ ~ 9co = 0.040. However, on increasing 
the lattice size the peak is shifted toward lower 9co values and the distribution function 
is clearly asymmetric. While the left-hand side remains practically independent of finite 
size effects (say for L > 128 LU), an increasingly long tail emerges on the right hand side 
of the distribution (for a detailed view see the inset of figure 8). The existence of these 
non- vanishing tails implies that fluctuations of the local CO coverage up to relatively 
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Figure 9: Plots of P^q versus L ^, where L is measured in LU. The straight hne corre- 
sponds to the best fit of the data that extrapolates to Pcoi^ ~^ °°) — 0.3544(2). 

large values are present in larger samples. Such excursions of the coverage induce the 
nucleation of critical CO clusters in larger lattices for smaller adsorption probabilities of 
CO, consequently Pqq{L) must decrease upon increasing L. In fact, figure 9, which shows 
a plot of Pqq{L) versus L^^, confirms this trend. As a first approximation the data can be 
fitted by a straight line that gives Pcoi^ ~^ — 0.3544(2). Furthermore our estimate 
for the coverage is 6qq = 0.043(1). These results point out that in the thermodynamic 
limit the spinodal point is very close to coexistence, i.e. AP^o = Pco ~ ^co^^ — 0.003. 
For the sake of comparison it is worth mentioning that for the ZGB model one has 



In order to perform epidemic studies it is necessary to account for the fact that the 
poisoned (absorbing) state above coexistence is non-unique, since it is due to a mixture 
of CO and N atoms with coverage 9co ~ 0.9 and On ~ 0.1, as shown in figure 1. 
So, the starting configuration has to be obtained running the actual dynamics of the 
system slightly above coexistence until 'natural' absorbing states suitable for the studies 
are generated. 

Extensive epidemic simulations ( ~ 2x10^ different runs) have been performed for the 
following values of Pco- Pco"'. Pco^ Pco{L = 256LU) and P^^(L = 256LC/). A value 
close to coexistence but slightly inside the active region, namely Pco = 0.347, has also 
been used. The obtained results can be observed in figure 10, which shows log-log plots 



AP 



CO = 



0.0012 [11]. 



3.2 Epidemic Study. 
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Figure 10: Log-log plots of the number of vacant sites N{t) versus t; measured in MCS; 
for epidemic simulations performed using lattices of size L — 256 LU. Results averaged 
up to 3x10^ different runs (triangle down) P^^ = 0.3544, (□) P§§ = 0.3522, {A)P§^^^ = 
0.35145, (O) Pco = 0-3506 and (0) Pco = 0.3470) . 



of N{t) versus t. It becomes evident that the method is quite sensitive to tiny changes 
of Pco- The obtained curves are fitted by equation (5), as shown in figure 10 for Pqo^ ■ 
The best fits are obtained for the parameters hsted in Table I. 

It is concluded that close to coexistence one has r^g// — 3.5 ± 0.5 with a characteristic 
time T ~ 50 MCS except for Pqo = 0.3506 with T ~ 160 MCS, a result that may reflect 
the fact that this point actually lies within the active region of the phase diagram. In 
fact, in the active region larger T values are expected, as judged by the results shown in 
figure 10. 

Of course, one may obtain better fits using different scaling ansatz than that in equa- 
tion (5). Our main finding, however, is that the occurrence of a power-law scaling be- 
havior close to coexistence can unambiguously be ruled out. This result is in qualitative 
agreement with recent numerical data obtained with the ZGB model [11]. All these ob- 
servations are also in agreement with the experience gained studying first-order reversible 
phase transitions where it is well established that correlations are short ranged, preventing 
the emergency of scale invariance. 

Let us now point our attention to the run performed taking Pqo = 0.347 (figure 10). 
A rapid decrease in N{t) up to 10^ MCS due to the low survivability of the initial empty 
patch is observed. However, for t > 10^ MCS only few epidemics survive (actually around 
10 epidemics over 10^) and the average number of empty sites increases according to 
N{t) oc t^, indicating the homogeneous propagation of the epidemic. For t ~ 10^ MCS 
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the active region covers the whole sample and the average density of empty sites remains 
stationary close to 9emp ~ 0.52. 



Pco 


Veff 


T MCS 


Number of epidemics 


= 0.3506 


4.2 ±0.2 


159 ±10 


4.3x10^^ 


p(^'oe^ = 0.35145 


3.45 ± 0.05 


58 ±1 


2.5x10^^ 


P^^ = 0.3522 


3.38 ±0.01 


50 ±1 


3.5x10^ 


= 0.3544 


4.1 ±0.1 


40 ±4 


5.5x10^ 



Table I Caption. Results obtained fitting the data corresponding to different values 
of Pco, shown in figure 10, using equation (5). 

4 Conclusions and outlook. 

In summary, after an extensive Monte Carlo simulation study of the first-order irreversible 
phase transition of the YK model for the catalized NO+CO reaction we concluded that: i) 
hysteretic effects arc found around the coexistence using the constant-coverage ensemble. 
The hysteresis loop basically exhibits three characteristic regions: a growing branch (GB), 
a vertical region (VR) and a decreasing branch (DB). Within each region the massive CO 
cluster corresponding to the inactive phase has a well defined convex (GB), concave (DB) 
and flat (VR) curvature. We are planning to study the dynamics of such interface in 
order to clarify the interplay between hysteresis and interfacial properties such as effective 
surface tension, roughness, curvature, etc. ii) The VR of the loop can be identified as 
the coexistence point, iii) The CC loop also gives evidence of a lattice size dependent 
upper spinodal point that can be extrapolated to the thermodynamic limit. Recently [20] 
one of us has shown that the spinodal point can be located quite accurately performing 
studies of the short-time dynamics of the system. So, we are planning to perform similar 
studies using the YK model in order to obtain an independent measurement of the upper 
spinodal point, iv) Epidemic studies reveal the existence of short ranged correlations close 
to coexistence. 

Based on these findings and the experience gained studying a similar system (namely 
the ZGB model [11]), we concluded that first-order irreversible phase transitions share 
many characteristics with their equilibrium (reversible) counterpart. So, we expect that 
these results will contribute to the development of a theoretical frame for the description 
of irreversible critical behavior. 
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